Rods of Neutron Scattering Intensity in Yb 2 Ti 2 07: Evincing a Hamiltonian with 
Significant Anisotropic Exchange in a Magnetic Pyrochlore Oxide 
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. Paramagnetic correlations in the magnetic material Yb2Ti207 have been investigated via neutron 

scattering, revealing a [111] rod of scattering intensity. Assuming interactions between the Yb 3+ ions 
composed of all symmetry-allowed nearest neighbor exchange interactions and long-range dipolar 
interactions, we construct a model Hamiltonian that allows for an excellent description of the neu- 
tron scattering data. Our results provide compelling evidence for significant anisotropic exchange 
interactions in an insulating magnetic pyrochlore oxide. We also compute the real space correlations 
leading to the [111] rod of scattering. 



0} ' In geometrically frustrated magnetic materials there exists no configuration of magnetic moments that simulta- 
^ ■ neously satisfies all the pairwise magnetic interactions. Experimental and theoretical research over the past twenty 
' years has shown that frustrated magnetic systems are prone to exhibit novel and intriguing collective thermodynamic 
phenomena 

Among frustrated three dimensional systems, the ^2^07 pyrochlores have attracted much attention 0. In these 
compounds, A is a trivalent rare earth ion (Ho, Dy, Tb, Gd, Yb) or yttrium (Y) and B is a tetravalent transition metal 
ion (Ti, Sn, Mo, Mn). Both A and B reside on two distinct lattices of corner-sharing tetrahedra. Theory predicts that 
classical Q and quantum Q Heisenberg spins on a pyrochlore lattice interacting via an isotropic antiferromagnetic 
nearest neighbor exchange Hamiltonian, Hn, fail to develop conventional LRO down to zero temperature. In real 
O . pyrochlore compounds, however, there generally exists some combination of other perturbing magnetic interactions 
(e.g. single ion anisotropy, dipolar interactions, etc) beyond Hr. Since Hn alone does not produce LRO, the low 
£SJ . temperature magnetic correlations of these materials are strongly influenced by the competition between material- 
specific perturbations. This is the origin of the richness of phenomena observed in the A2-B2O7 pyrochlores 
including spin liquid spin glass Q, spin ice Q, and LRO with persistent low-temperature spin dynamics 
In this article, we consider the Yb2Ti207 pyrochlore which does not apparently exhibit any of the aforementioned 
1^ . phenomena and has some unique and unusual features of its own which have heretofore remained unex pla ined. 

Yb2Ti207 has a ferromagnetic character with a Curie-Weiss temperature, 6*cw = +0.65 ± 0.15 K [IE ED]- The 
Yb 3+ ~ 3 /ib magnetic moments predominantly lie perpendicular to the local [111] cubic unit cell diagonals, making 
this system the only known local [111] XY pyrochlore with a ferromagnetic Ocw @]- Magnetic specific heat (C m ) 
measurements reveal a sharp first order transition at T c ps 240 mK [l^, suggesting the onset of LRO. While a single 
crystal elastic neutron scattering (NS) study suggested ferromagnetic order below T c [3], a subsequent polarized 



NS study [L| did not confirm such ordering. Furthermore, powder NS shows no LRO down to 110 mK [151 ] and 
very recent NS on a single crystal sample has not found any sign of LRO in a broad region of the (hkk) scattering 
plane at 30 mK [l6jj. The T c « 240 mK transition seen in C m has therefore so far not been matched with the 
observation of conventional (dipolar magnetic) LRO. In addition, Mossbauer spectroscopy and muon spin relaxation 
(/iSR) measurements find a rapid decrease of the Yb 3+ magnetic moments fluctuation rate, v, upon approaching T c 
from above, with //SR revealing a temperature-independent v (i.e. persistent spin dynamics) from T c down to 40 
mK, the lowest temperature considered [l5[. Considering all these results together, one may ask whether the 240 
mK transition in Yb2Ti207 may be another rare example of hidden (non-dipolar) order [171 ]. Another intriguing 
possibility [HI is that the 240 mK first order transition takes place between a "spin gas" (paramagnetic) state and a 
spin liquid without any symmetry breaking. 

A very interesting feature of the magnetic correlations in Yb2Ti207 found at temperatures T c < T < 2 K are rods 
of NS intensity along the [111] directions 0, [l6|. At first sight, the presence of such rods signals an anisotropy in 
the magnetic correlations that may originate from a structural transition at T > 2 K or from intrinsically stronger 
correlations within the kagomc planes perpendicular to the four [111] directions forming the undistorted pyrochlore 
structure as compared to correlations perpendicular to the kagome planes [l6| — hence making a "spin liquid crystal" 



2 



of sorts. 

Here, we report results from diffuse NS measurements on Yb2Ti207 in a temperature range above 9cw > T c . By 
numerically annealing a set of exchange couplings to maximize the agreement between experimental results and NS 
computed within a random phase approximation (RPA), we determine a spin Hamiltonian, H, that captures the 
main features of the observed paramagnetic NS pattern and reveals that significant spin exchange anisotropy exists 
in this insulating pyrochlore oxide material. There have been recent claims of evidence for anisotropic exchange 
at play in Yb2Ti2C>7 [HI, [l9[ and other A2B2O7 pyrochlores 18, 2(J. However, because of the limitations imposed 



by the physical quantity considered (local susceptibility), and the models used in these works, the specific nature 



and symmetry of the putative microscopic exchange has, until this work remained hidden [21|. Owing to the highly 
structured spin correlation functions in Yb2Ti207, we obtain compelling evidence for anisotropic exchange. From the 
reciprocal space structure factor, we compute the spin-spin correlations along different crystallographic directions. 
While these correlations are anisotropic, the correlation lengths themselves do not distinguish between correlations 
parallel and perpendicular to the kagomc planes. 

The NS cross-section was measured on the D23 diffractometer at the Institut Laue Langevin, France. A single 
crystal rod was aligned with [hOO] and [Okk] in the scattering plane, hence providing access to all principal symmetry 
directions of the cubic crystal structure. With incident neutron energy of 14.7 meV, significantly larger than any 
characteristic energy scale in the system, and counting all final energies, the measured intensity is proportional to the 
spatial Fourier transform of the instantaneous correlation function 5(q) = J S(c[,uj)duj 

Figs. CE{a,b) show experimental NS data for Yb 2 Ti 2 7 at T = 9.1 K and T = 1.4 K in the (hkk) plane. Sharp 
intense Bragg peaks at integer reciprocal lattice positions were removed from the data to expose significant structure 
in the diffuse magnetic scattering. Fig. QJi shows the NS map at T = 9.1 K. The magnetic correlations weaken with 
increasing temperature. Indeed, most of the features present at T = 1.4 K (Fig. [TJd) arc absent at 9.1 K, with only 
a weakened rod of scattering along [111] and a feature in the upper right corner at 3.5, 2.25, 2.25 remaining. Fig. [TJd 
shows the NS map at T = 1.4 K, where the most interesting feature is the aforementioned rod of scattering along 
the [111] direction 0,[lB]. Figure [Dd exhibits other features of interest such as a weaker rod of scattering going from 
400 to 222 and intensity near the point 022. The intensity of the feature centered on 3.5, 2.25, 2.25 does not change 
with temperature (Fig. QJ) indicating that it is not magnetic in origin and can therefore be omitted from further 
consideration. 

To explain the NS pattern, we propose a Hamiltonian, H = H c { + Hi nt , that includes a crystal field (CF) part, H c { , 
and spin-spin interactions, Hi nt = H^ip + H cx . The form of H c { is fixed by the symmetry of the Yb 3+ environment. 
The two sets of CF parameters that we use have been determined in Refs. [la. Hjjj]. The magnetic Yb 3+ ion has 
electronic configuration 2 F 7 / 2 , hence J = 7/2 and Lande factor g,j =8/7. The nearest neighbor distance between 



Yb 3+ ions is r, m = (a/4)v / 2, where a = 10.026 A is the size of the conventional cubit unit cell 14|. This fixes the 



strength of the coupling D = ^j^r^ra ~ 0.01848 K of the long-range magnetostatic dipolar interaction, 
5 Efi a-,i b) D m?$ ( J ? ' J ? - 3 W ' R-tf )( J i ' R« ))■ Wc also consider H ox which contains all nearest neighbour exchange 
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interactions, J ei that respect lattice symmetries. There are four such nearest neighbor interactions [22|: -ffismg = 
- Jising S<i, y,6> ( J ? ' ( J ? ' 2& ) ' which couples the local [111] z components of J, H iso = -J iso V , „ :; . J® • J^, 
the standard isotropic exchange, H pd = —J p & J2<i a -j b>(^i ' ^ ~ ' ' -^y ))i a pseudo-dipolar interaction 

of exchange origin and not part of H d i P and, finally, #dm = -Jdm J2<i a ;j b> ^dm ' (^i x tne Dzyaloshinskii- 
Moriya (DM) interaction [23|. In all of these terms, Jf denotes the angular momentum of the Yb 3+ located at lattice 
R° (FCC lattice site i, and tetrahedral sub-lattice site a) [25| and R^ 6 is a unit vector directed along — R". The 
relationship between H- lnt and a corresponding effective spin- 1/2 model is discussed in Ref. (24| . 

We use H to compute the diffuse NS pattern within the RPA 2(| 27] (see Ref. 24] for justification of the usage of 



RPA here). We first compute the single ion susceptibility, x^°\ from H, 
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where n v is the thermal occupation fraction for CF state v. M" = Eg( t/ |J a lA*) M a,aj where u% a is the rotation 
matrix from the local (a) frame defined on sublattice a to the global (a) frame. The operator 3 a acts on the CF states 
defined in the local a quantization frame. The CF wavefunctions |z/), at sublattice site a, are obtained by diagonalizing 
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FIG. 1: (Color online) Neutron scattering (NS) maps in the (h,k,k) plane, (a) and (b) show experimental data at 9.1 K and 
1.4 K, respectively. In (b), a rod of scattering intensity along [111] is clearly seen, while in (a) the rod is not strong, (c) shows 
a comparison between experimental data and the computed NS intensity at 1.4 K along [3 — k, k, k] for both isotropic exchange 
and the full model, (d) shows computed NS using the Hamiltonian H (see text), at T = 1.4 K. (e) shows the calculated NS 
using H with only isotropic exchange determined from fitting #cw and long-range dipolar interactions at T = 1.4 K. (f) shows 
experimental neutron scattering at 1.4 K and 9.1 K along the line [h, 2.25, 2.25] showing that the feature between h = 3 and h = 
4 does not change in intensity with temperature. The white arrows in (a) and (b) indicate the range h £ [3, 4] in [h, 2.25, 2.25] 
in (f). 

ffcf 00. The interacting RPA susceptibility, x (q 1 w) ™ then (q, u) + E 7Ac X°' a ^)J f(q)x^ = 
5 ab Xa' a ^ {u), where J(c\) is the Fourier transformation of the interaction matrix J{i,j) for Hamiltonian iJi nt = 
— (1/2) J2i j- a b a p ^ta^ab (^j)^j b- The infinite lattice sum of the dipolar interaction is computed using Ewald sum- 
mation (25[. We solve for (q, ^) numerically. Finally, the NS function, 5(q,w) (26, 27|, is given by 
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where / (Q) is the magnetic form factor for the Yb 3+ ion [231. Q = q + G is the scattering wave vector where q is a 
wave vector inside the first Brillouin zone, and G is an FCC reciprocal lattice vector. r a and r b are basis vectors for 
the tetrahedral sublattice (25j |. 

The fit to the experimental data was performed by computing the RPA scattering intensity along the measured lines 
in Q space (dashed lines in Fig. [T|b). The computed intensities were rescaled to the experimental count rate using 
the relation S' (q) = c.qS (q) + c\ + C2 | Q | [3Cj| . where the parameters Co, c\, and C2 are the same for all Q points. 
There arc therefore seven adjustable parameters in total with the four exchange couplings, J e and the three c„ fitting 
parameters. The variance between the measured and calculated neutron data, with a contribution to the variance 
from fitting #cw as well, was minimized using a simulated annealing algorithm |3lj |. As a first step, for simplicity 
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and computational speed, we make a static approximation [32| to 5 (q, w). Within the moderate constraints 24 of 
the procedure followed, the present calculation is suitable to reach the main conclusion of this work: that significant 
anisotropic exchange couplings J £ are necessary to account for the structure of the NS pattern of Yb 2 Ti 2 07. 

Fig. [TJl shows the RPA NS pattern in the (hkk) plane at 1.4 K obtained from simulated annealing fits to the 
experimental neutron scattering map of Fig. [Tp. The model data in Fig. [T}i is obtained from a H ex with J7i s i n g = 0.81 
K, Jiso = 0.22 K, J7pd = —0.29 K, and Jdm = —0.27 K. The calculated intensity matches the experimental data 
(Fig. [Tp) well, providing strong evidence that our model H cx contains the correct interactions for Yb2Ti2C>7. A 
cut along [3 — k, k, k] (Fig. [TJ:) emphasizes the quantitative agreement between the computed and experimental NS 
intensities. A similar quality of fit was obtained for other line scans (dashed lines in Fig.QJ>) shown in the supplemental 
material [13]. We carried out the fitting procedure using CF parameters for H c { taken from Refs. 15, 3] finding that 
the exchange couplings do not change significantly. Fig.[Tp and the dashed (black) line in Fig.[Tp show the NS intensity 
calculated for a model H lrA with only long-range dipolar and isotropic exchange interactions (Rising = Jpd = Jr>u = 0) 
with J- lso = 0.06 K, determined by fitting #cw, and refitted c„'s. Clearly this model and the resulting NS pattern 
do not describe the experimental data (Fig. [TJ)) well at all. Similarly, a Hamiltonian with isotropic-only exchange 
does not describe the local susceptibility, xiocai, well lB"21|. On the other hand, the present anisotropic Hamiltonian 
describes xi OC ai with no adjustable parameters [2~j| . 

To rationalize the direct space origin of rods of NS intensity, we computed the spin-spin correlation function from the 
reciprocal space RPA susceptibility (J" (r) (0)) = k^,T J x"b ex P ' r ) where r = R^ — Rf . The integral was 
performed numerically over the first Brillouin zone. We considered the isotropic real space correlations (by summing 

over all directions a in spin space) S(r) = ^2 a (Ja ( r ) Jj? (0)), and also S±(r) = k^T f[S a p — 4a4p]Xab (l) CX P (*1 ' r ) ^ 
whose Fourier transform is measured in the NS. Figure [5] shows S(r) and S±(r) for r taken along the [111] direction, 
and [011], [121] perpendicular to [111]. S(r) and S±(r) are similar, so we discuss both together. The correlation 
lengths for these three directions were extracted, assuming exponential decay, and found not to differ greatly within 
our margin of error. However, the [011] correlations are larger than those of the other two directions indicating 
some degree of spontaneous decoupling of the kagome planes albeit perhaps with a quasi-isotropic correlation length. 
Since the [Oil] directions lie within two kagome planes and the [121] directions lie in only one, truly two dimensional 
correlations may have been be expected to lead to the splitting S(r)r i]i > S , (r)[i2i] > ^(^[xii]- This is not borne out 
by our results — the similarity of S(r) for the [121] and [111] directions indicates that the correlations are strongest 
along spin chains. Whereas the NS intensity in the (hhk) plane, considered on its own, suggests that the correlations 
are quasi two-dimensional flij |. with a weak decoupling of the kagome planes, the real space correlations within our 
model do not support this simple picture. 

From the determined H ox , RPA predicts a second order phase transition to a ferromagnetic phase (ordering wavevec- 
tor q = 0) at a critical temperature T C RPA w 1.2 K. We expect thermal and quantum fluctuations to renormalize the 
values of the anisotropic exchange J e determined above. 




FIG. 2: (Color online) Real space spin-spin correlation functions (a) S(r) and (b) S±(r) (see text) computed at T = 1.4 K 
using the anisotropic exchange model and plotted along various crystallographic directions. 

In summary, we have presented diffuse neutron scattering (NS) maps of Yb2Ti2C>7 in its paramagnetic regime, 
finding rods of scattering in the (111) directions. By fitting this data to NS computed from a candidate model Hamil- 
tonian, we found a set of couplings that reproduce the main features of this NS pattern. The Hamiltonian includes, 
as an essential component, sizeable anisotropic exchange interactions. This suggests that anisotropic exchange might 
be important in other .A2.B2O7 rare earth magnets 3 2(|. We find that the rods of scattering occur without any 
symmetry breaking from, for example, a structural phase transition. We anticipate that our results will allow for a 
greater understanding of the nature of the phase transition at 240 mK and of the low temperature phase of Yb2Ti207. 

We thank Y.-J. Kao for useful discussions and contributions at the earliest stage of this project. We acknowledge 
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SUPPLEMENTAL MATERIAL 



This document is a supplement to our main article, where we explain the relationship between our model Hamilto- 
nian and effective spin-1/2 models appropriate for this material, such as the model considered in Ref. [33| for P1-2A/2O7 
(M = Ir, Zr, Sn) pyrochlore materials. We also provide details on the validity of the random phase approximation 
(RPA) for computing the neutron scattering pattern at T = 1.4 K, constraints on the results of our RPA calculations, 
a short discussion of energy integrated neutron scattering results, and additional line scans to demonstrate the success 
of our model. 



EFFECTIVE SPIN HALF EQUIVALENT HAMILTONIAN 

In Yb2Ti2C>7, the energy gap separating the crystal field ground state and the first excited doublet is A ~ 620 K, 
very large compared to the magnetic dipolar and exchange (bilinear or higher multipolar |17| ) interactions at play 



in the compound. Hence, as in spin ices [34| and LiHoF 4 [35|, but unlike the Tb2Ti 2 07 [28§ and Tb 2 Sn 2 07 [36[ 
pyrochlores, it is safe to project the microscopic Hamiltonian, H lnt , describing the interactions between Yb 3+ ions, 
into the Hilbert space spanned solely by a direct product of the states in the Yb 3+ crystal field ground state doublet. 
Such a projection allows us to rewrite H{ nt , irrespective of its original complicated (even possibly multipolar) form, in 
terms of an effective Hamiltonian, -ff e ff, expressed solely in terms of anisotropic couplings between effective spin-1/2 
operators. In principle, we could have considered such an effective model in the main body of our paper, parametrized 
by effective couplings jTiso, Rising, J7pd and J7dm, and expressed the neutron scattering function in terms of components 
of J projected onto the crystal field ground state doublet. Again, this is possible since there is negligible operator 
correction to the two-point correlation function that enters in the scattering function because A is so large compared 
to H int . 

In this context, it is therefore important to view the model describing the exchange interactions, H cx , with couplings 
Jiso, Rising, Jpd and Jum, not as a microscopic model capturing the correct physics on energy scales comparable to 
the crystal field splitting (which may include multipolar interactions) [l7| . but rather one whose projection gives the 
correct effective low-energy theory in terms of pseudo-spin- 1/2. Our model is therefore an "un-projected" version of the 
effective spin-1/2 Hamiltonian H c jf describing Yb2Ti207, which we have taken to be a strictly bilinear Hamiltonian, 
^F£j v JfJJ. This was done for presentation sake and to relate our approach to what was used previously in a random 



phase approximation treatment of Tb2Ti207 27j. In the rest of this section, we relate our bilinear couplings {J7e} to 
the anisotropic couplings between the components of a spin-1/2 in an effective i/ e ff Hamiltonian. 

The exchange model H cx in our work consists of bilinear exchange terms between the full blown angular momentum 
operators J for the Yb 3+ ions. As stated above, the exchange model can be written in terms of an effective spin-1/2 
model. This is done by computing the expectation value of the three components of J within the ground state doublet 
and defining g\\ = 2gj(tp± \ J z \ i/j±) and g± = 2gj(ip± | J x \ ijj±) = 2g.j(ip± | 3 V | tp±)- Here | ip±) are the two states 
that make up the ground state doublet of the crystal field, J a is expressed in the local coordinate system for each 
corner of a tetrahedron, gj = 8/7 is the Lande factor, and z is the corresponding [111] cubic diagonal direction. In 
this coordinate system, g\\ and g± are the only non-zero terms, and we can construct the g tensor: 

(3) 

For the crystal field of Ref. [ll( (Ref. [l9|), g\\ = 1.77 {g\\ = 2.25) and g ± = 4.18 (g x = 4.1). Using this g tensor, 
the relationship between the full blown angular momentum operators and the effective spin-1/2 operators is given by 
J = ■fjSeff = 29?^' wnere "Seff is an effective spin-1/2 and a is a vector of Pauli matrices. 

Using this relationship, we can recast H ex in terms of effective spin-1/2 operators by replacing all of the terms Jf 

with j£j<?i where the labels i and a are required because we arc now passing from a local coordinate system to the 
global cartesian coordinate system, meaning that g" now changes from site to site. Performing this transformation, 
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we obtain H cS = H[ sing + H[ so + H' pd + H^ M , where: 

flU = -^ £ E way-**) (4) 

^ <i,a-,j,b> 

H L = -%- E bW-sM. (6) 

J <i,a:j,b> 

H U = -%* E (gK-gM-«'^)(g?4-^)). ( 6 ) 

<i,a;j,b> 

4m = ~^ E "dm ■ (g?a? x §M) • (7) 

J <i,a;j,b> 

In addition to the representation of the symmetry allowed bilinear exchange interactions used in this paper, other 
representations are possible, which are linear combinations of the terms used in the original representation of H e g. 
One other representation is given by Ref . [37} : 



(r,r') 



nn 

Jim E \^KK' +9 X iplK> +&rK>) 



+g q y { &r ■ Tt r ,r' j { °V' ■ Tir,r' ) — ( &r ' n r . 
+g K I I a r i ■ n r _ r > I + I a T ■ n r _ r > \ d* T 



(8) 



In this representation x, y, and z refer to the local coordinates at each corner of the tetrahcdra, oy = (<T^,cr^), 
n r ,r' = (cos^r,r'5 — sin0 r r f), rij, r , = (sin</> r ,r') cos</v,r')j and </> r . r ' = 0, 27r/3, — 27r/3 (37| |. The various g a terms arc 
defined as: 

S« = 1 - sV&r - |t 2 - zVGx 3 + ^x 4 , (9) 

.9^ = 1 + iVEx + ^-x 2 - 3V6x 3 + ^-x 4 , (10) 
2 16 

g q =-2(l- 2v / 6a- + 9a; 2 - 3\/6a; 3 + ^x 4 ^ , (11) 

g K = 2V2 (l + V6x - ^z 2 + - ^x 4 ) , (12) 

where x = Vpf^/Vpfa, is the ratio of two Slater-Koster parameters, representing transfer integrals between p x /p y and 
fx(5z 2 -r 2 )/ fy(5z 2 -r 2 ) orbitals and p z and /(5 Z 2_ 3r 2) z orbitals, respectively [37}. This representation can be related to 
our original notation by using the relationships between global cartesian coordinates and the local coordinate system 
for each corner of the tetrahedral sublattice, and the relationships between the exchange terms in the body of the 
paper and the exchange terms in Eq. [5] The net result of these relationships is expressed by the following relations: 

- Jnn.9 1 ' = MJlsing + Jiso - 5Jpd - 4J DM ) , (13) 

-Jnng ± = ~^jf^2 (^ is ° _ 2^ pd + 2^ DM / ' 

9l "L , 7 



Jnng 9 — 24^2" [^iso + J^pd — JdM ) , (15) 



'Jm\g K — \2g^ — (^Jiso — 2^P d ^»7dM^ 



(16) 



By inverting these relations, we can compute the ratios ^ i3 ° , and j° M , from the work of Ref. 37 1, finding 

that, for x w [—00, —0.5] and a; « [20, 00], there is fair agreement in terms of sign and magnitude with those of our 



■Jmo vTpd v7dM 

Rising Rising Rising 

FIG. 3: (Color online). Plots of the ratios j^' B ° , J^ pd , and ^ DM as a function of x computed from the work of Ref. [3?| for 
the ranges x « [—20, —0.5] and x « [10, 30]. 



work = 0.28, = -0.36, and #M. = -0.33, sec Fig. O In the range x w [-0.5, 101, the ratios 

Rising ' Rising Rising O 1— f O L ' J ' Rising' JDM ' 

and ^ DM diverge when Jis^g — > 0,and fluctuate rapidly, making any comparison to our model very difficult. 

In addition to these effective spin-1/2 models, we can relate the notation used in this work, H ex = £zi s i ng + H lso + 
H p d + ffDM, to the form of the nearest neighbour invariants described in Ref. [22| as {%}• These are: 

xi = + Jp y 2 + Jir 3 + jfjg + j? jf + (17) 
+J 2 j| + J2J3 + J2 J4 + J2J4 + J3 J4 + J3 J4; 

X2 = JfJ|+JpI + J?^ + J^ + JK + JfJI, (18) 
X3 = JfJl + J?Jf + J?Jl+JfJ| + J?JI + JfJl (19) 

T»/ TZ TZ TV _ Till _ TZ TI _ TI t!/ _ fH Tl 

d 2 d 3 d 2 d 3 J 2 J 4 d 2 d 4 J 3 J 4 d 3 J 4> 

X 4 = [J?, JI] + [J? , JI] + [J?, Jg] + [Jf , Jg] + Pi, Jf ] + [J?, J2] (20) 
+ [Jf , Jf] + [jf , Jf] + [Jf , jf] + m , JI] + [Jf, Jf] + [Jf , Ji] , 

where [J", Jj] = J" JJ — Ji Jf - In this notation, we may write -Hex = — — ^2X2 — >/3X3 — JiXi- The relationships 
between the {Je\ and { J n } are given by: 

Jlsing = -3J!+3J2+3J 3 (21) 

Jiso = ^ J l + ^ J 2 + j J 3 (22) 

2 2 4 

Jpd = - g J2 - 2^3 (23) 
Jdm = J1-J2-J3 + J4 (24) 

JUSTIFICATION FOR THE USE OF THE RANDOM PHASE APPROXIMATION 

The random phase approximation (RPA), which uses the bare single- ion anisotropy as non- interacting reference 
susceptibility, is well justified in the case of Yb2Ti 2 07 at T = 1.4 K for two main reasons. Firstly, H ex and ifdip 



are weak (~ 10° K and ~ 0.02 K, respectively) compared to the large (A ~ 620 K llj) energy gap between the 
ground state crystal field doublet and the first excited crystal field doublet [HI]. This implies that there is negligible 
interaction-induced admixing between the ground and first excited crystal field doublets, unlike the case of Tb2Ti2C>7, 
where admixing between the two lowest energy crystal field doublets is significant, due to the much smaller energy gap 
between between the ground and first excited states (~ 18 K) [2^|. Secondly, the correlations used to calculate the 
neutron scattering are at 1.4 K, a temperature of about twice the Curie- Weiss temperature (0cw ~ 0.75 K 13, HI]), 
with the system in the paramagnetic regime, where RPA should be reasonably valid. 
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CONSTRAINTS OF THE RPA METHOD AND STATIC APPROXIMATION 

For the set of parameters, J e = {Rising, J\so, J-pd, iTdm}, obtained from our minimization procedure using the static 
approximation, we carried out a more rigorous energy integrated calculation of the full dynamical scattering function, 
S (q, uj). Performing this calculation, we obtain a very similar, but not quantitatively identical, reciprocal space map 
of intensities to the scattering pattern obtained using the static approximation. 

When using the RPA, it is also important to consider whether or not correlation effects, that develop near phase 
transitions, are being neglected. In this work we use the RPA to perform calculations at T = 1.4 K, approximately 
5 times the reported transition temperature of Yb2Ti2C>7 (~ 240 mK (HI), an d twice the Curie- Weiss temperature 
(6*cw ~ 0.75 K [13, EH)- Due to these correlation effects, we expect the fitted anisotropic exchange coupling to 
be renormalized as a function of the temperature at which the fit is computed. The change should be of order 
1/ (T — #cw) 2 ~ 25%, the leading correction in a high temperature expansion of x(q), compared to that of the 
RPA approximation of \ (q)- Despite these constraints, each unique set of couplings will generate a distinct diffuse 
scattering pattern that should be captured well by the RPA calculations to within a scaling factor that will differ 
from the scaling factor determined from the RPA calculations by a factor of order one. This is based on the fact that 
the physical scale is fixed by the Curie- Weiss temperature, and we included this as a constraint when performing the 
simulated annealing fit. 

ADDITIONAL NEUTRON SCATTERING LINE SCANS 

To supplement the results presented in Fig. Q] of the main body of the paper, Fig. 2] shows additional cuts through 
the [hkk] plane at T = 1.4 K. We chose this particular temperature in an effort to maximize the signal to noise 
ratio of the [111] rod feature. Using data collected at higher temperatures, where the RPA would suffer less from 
temperature renormalization effects, would mean a loss of scattering intensity, and thus greater difficulty distinguishing 
the features of interest from background scattering. Higher temperature data with high signal-to-noise ration would 
be very desirable to ascertain further the quantitative accuracy and aforementioned temperature renormalization of 
the exchange parameters determined in the present work. 
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FIG. 4: (Color online). Panels a through f show linescans within the [hkk] plane. The blue dots are experimental diffuse neutron 
scattering data taken at T = 1.4 K. The solid red lines are calculated from the single set of anisotropic exchange couplings 
( Jiso, Jising, Jpd, Jom) presented in the body of the paper, and the black dashed line is the prediction of a model with isotropic 
exchange and dipolar interactions, with the scale of the isotropic exchange set solely by the value of the Curie- Weiss temperature 
(6cw ~ 0.75 K [13, El) 



10 



H.T. Diep, Frustrated Spin Systems, (World Scientific Publishing Co. Pte. Ltd., 2004, Hackensack). 

J. S. Gardner, M. J. P. Gingras, and J. E. Greedan, Rev. Mod. Phys. 82, 53 (2010). 

R. Moessner and J. T. Chalker, Phys. Rev. Lett. 80, 2929 (1998). 

B. Canals and C. Lacroix, Phys. Rev. Lett. 80, 2933 (1998). 

J. S. Gardner et al, Phys. Rev. Lett. 82, 1012 (1999). 

M. J. P. Gingras et al, Phys. Rev. Lett. 78, 947 (1997). 

S. T. Bramwell and M. J. P. Gingras, Science 294, 1495 (2001). 

Y. Chapuis et al., Physica B 404 686 (2009). 

P. Bonville et al., Hyperfine Interactions 156, 103 (2004); ibib. |arXiv:cond-mat/0306470| 
S. T. Bramwell et al, J. Phys.: Condens. Matter 12 483 (2000); 
J. A. Hodges et al., J. Phys. Condens. Matter 13, 9301 (2001). 

H. W. J. Blote, R. F. Wielinga, and W. J. Huiskamp, Physica (Amsterdam) 43, 549 (1969). 

Y. Yasui et al, J. Phys. Soc. Jpn. 72, 3014 (2003). 

J. S. Gardner et al., Phys. Rev. B 70, 180404(R) (2004). 

J. A. Hodges et al, Phys. Rev. Lett. 88, 077204 (2002). 

K. A. Ross et al., Phys. Rev. Lett. 103, 227202 (2009). 

P. Santini et al, Rev. Mod. Phys. 81, 807 (2009). 

H. Cao et al, Phys. Rev. Lett. 103, 056402 (2009). 

H. B. Cao et al., J. Phys.: Condens. Matter 21, 492202 (2009). 

B. Z. Malkin et al., J. Phys.: Condens. Matter 22, 276003 (2010). 

J. D. Thompson et al. J. Phys.: Condens. Matter 23, 164219 (2011). 

P. A. McClarty et. al., J. Phys.: Conference Series 145, 012032 (2009). 

M. Elhajal et al. Phys. Rev. B. 71, 094420 (2005). 

See the Supplemental Material. 

M. Enjalran and M. J. P. Gingras, Phys. Rev. B 70, 174426 (2004). 

J. Jensen and A. R. Mackintosh, Rare Earth Magnetism, (Clarendon Press, 1991, Oxf ord). 
Y. -J. Kao et al, Phys. Rev. B 68, 172407 (2003). 

H. R. Molavian, M. J. P. Gingras and B. Canals, Phys. Rev. Lett. 98, 157204 (2007). 

P. J. Brown in International Tables for Crystallography, edited by A. J. C. Wilson (Kluwer Academic Publishers, Dordrecht, 
1992), Vol C. p 393. 

T. Yavors'kii et al., Phys. Rev. Lett. 101, 037204 (2008). 

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery Numerical Recipes: The Art of Scientific Computing, 
(Cambridge University Press, 2007, Cambridge). 

S. W. Lovesey, Theory of Neutron Scattering from Condensed Matter (Clarendon Press, 1984, Oxford). 
S. Onoda and Y. Tanaka, Phys. Rev. B 83, 094411 (2011). 

M. J. P. Gingras, in Introduction to Frustrated Magnetism Materials, Experiments, Theory, edited by C. Lacroix et al. 
Springer Series in Solid-State Science Vol. 164 (Springer- Verlag, Berlin, 2011). 

S. M. A. Tabei et al. Phys. Rev. B 78, 184408 (2 008). 

P. A. Mc Clarty, P. Stasiak, and M. J. P. Gingras larXiv: 101 1.63461 (201 1) . 
S. Onoda, larXiv: 1101. 1230 1 



